function y = FUN_whitenoise(average,variance,size)

y = randn(1,size); 
y = y - mean(y); 
y = y/std(y); 
std_variance = sqrt(variance); 
y = average + std_variance * y;

end